Parameter and coupling estimation in small networks of Izhikevich's neurons.

Nowadays, experimental techniques allow scientists to have access to large amounts of data. In order to obtain reliable information from the complex systems that produce these data, appropriate analysis tools are needed. The Kalman filter is a frequently used technique to infer, assuming a model of the system, the parameters of the model from uncertain observations. A well-known implementation of the Kalman filter, the unscented Kalman filter (UKF), was recently shown to be able to infer the connectivity of a set of coupled chaotic oscillators. In this work, we test whether the UKF can also reconstruct the connectivity of small groups of coupled neurons when their links are either electrical or chemical synapses. In particular, we consider Izhikevich neurons and aim to infer which neurons influence each other, considering simulated spike trains as the experimental observations used by the UKF. First, we verify that the UKF can recover the parameters of a single neuron, even when the parameters vary in time. Second, we analyze small neural ensembles and demonstrate that the UKF allows inferring the connectivity between the neurons, even for heterogeneous, directed, and temporally evolving networks. Our results show that time-dependent parameter and coupling estimation is possible in this nonlinearly coupled system.

Inferring the connectivity of coupled chaotic oscillators using Kalman filtering.

Inferring the interactions between coupled oscillators is a significant open problem in complexity science, with multiple interdisciplinary applications. While the Kalman filter (KF) technique is a well-known tool, widely used for data assimilation and parameter estimation, to the best of our knowledge, it has not yet been used for inferring the connectivity of coupled chaotic oscillators. Here we demonstrate that KF allows reconstructing the interaction topology and the coupling strength of a network of mutually coupled Rössler-like chaotic oscillators. We show that the connectivity can be inferred by considering only the observed dynamics of a single variable of the three that define the phase space of each oscillator. We also show that both the coupling strength and the network architecture can be inferred even when the oscillators are close to synchronization. Simulation results are provided to show the effectiveness and applicability of the proposed method.

Consistency of heterogeneous synchronization patterns in complex weighted networks.

Synchronization within the dynamical nodes of a complex network is usually considered homogeneous through all the nodes. Here we show, in contrast, that subsets of interacting oscillators may synchronize in different ways within a single network. This diversity of synchronization patterns is promoted by increasing the heterogeneous distribution of coupling weights and/or asymmetries in small networks. We also analyze consistency, defined as the persistence of coexistent synchronization patterns regardless of the initial conditions. Our results show that complex weighted networks display richer consistency than regular networks, suggesting why certain functional network topologies are often constructed when experimental data are analyzed.

Integration of cellular signals in chattering environments.

Cells are constantly exposed to fluctuating environmental conditions. External signals are sensed, processed and integrated by cellular signal transduction networks, which translate input signals into specific cellular responses by means of biochemical reactions. These networks have a complex nature, and we are still far from having a complete characterization of the process through which they integrate information, specially given the noisy environment in which that information is embedded. Guided by the many instances of constructive influences of noise that have been reported in the physical sciences in the last decades, here we explore how multiple signals are integrated in an eukaryotic cell in the presence of background noise, or chatter. To that end, we use a Boolean model of a typical human signal transduction network. Despite its complexity, we find that the network is able to display simple patterns of signal integration. Furthermore, our computational analysis shows that these integration patterns depend on the levels of fluctuating background activity carried by other cell inputs. Taken together, our results indicate that signal integration is sensitive to environmental fluctuations, and that this background noise effectively determines the information integration capabilities of the cell.

Relating structural and functional anomalous connectivity in the aging brain via neural mass modeling.

The structural changes that arise as the brain ages influence its functionality. In many cases, the anatomical degradation simply leads to normal aging. In others, the neurodegeneration is large enough to cause neurological disorders (e.g. Alzheimer's disease). Structure and function can be both currently measured using noninvasive techniques, such as magnetic resonance imaging (MRI) and electroencephalography (EEG) respectively. However, a full theoretical scheme linking structural and functional degradation is still lacking. Here we present a neural mass model that aims to bridge both levels of description and that reproduces experimentally observed multichannel EEG recordings of alpha rhythm in young subjects, healthy elderly subjects, and patients with mild cognitive impairment. We focus our attention in the dominant frequency of the signals at different electrodes and in the correlation between specific electrode pairs, measured via the phase-lag index. Our model allows us to study the influence of different structural connectivity pathways, independently of each other, on the normal and aberrantly aging brain. In particular, we study in detail the effect of the thalamic input on specific cortical regions, the long-range connectivity between cortical regions, and the short-range coupling within the same cortical area. Once the influence of each type of connectivity is determined, we characterize the regions of parameter space compatible with the EEG recordings of the populations under study. Our results show that the different types of connectivity must be fine-tuned to maintain the brain in a healthy functioning state independently of its age and brain condition.

Relaxation dynamics and frequency response of a noisy cell signaling network.

We investigate the dynamics of cell signaling using an experimentally based Boolean model of the human fibroblast signal transduction network. We determine via systematic numerical simulations the relaxation dynamics of the network in response to a constant set of inputs, both in the absence and in the presence of environmental fluctuations. We then study the network's response to periodically modulated signals, uncovering different types of behaviors for different pairs of driven input and output nodes. The phenomena observed include low-pass, high-pass, and band-pass filtering of the input modulations, among other nontrivial responses, at frequencies around the relaxation frequency of the network. The results reveal that the dynamic response to the external modulation of biologically realistic signaling networks is versatile and robust to noise.

Nonlinear chemoconvection in the methylene-blue-glucose system: two-dimensional shallow layers.

Interfacial hydrodynamic instabilities arise in a range of chemical systems. One mechanism for instability is the occurrence of unstable density gradients due to the accumulation of reaction products. In this paper we conduct two-dimensional nonlinear numerical simulations for a member of this class of system: the methylene-blue-glucose reaction. The result of these reactions is the oxidation of glucose to a relatively, but marginally, dense product, gluconic acid, that accumulates at oxygen permeable interfaces, such as the surface open to the atmosphere. The reaction is catalyzed by methylene-blue. We show that simulations help to disassemble the mechanisms responsible for the onset of instability and evolution of patterns, and we demonstrate that some of the results are remarkably consistent with experiments. We probe the impact of the upper oxygen boundary condition, for fixed flux, fixed concentration, or mixed boundary conditions, and find significant qualitative differences in solution behavior; structures either attract or repel one another depending on the boundary condition imposed. We suggest that measurement of the form of the boundary condition is possible via observation of oxygen penetration, and improved product yields may be obtained via proper control of boundary conditions in an engineering setting. We also investigate the dependence on parameters such as the Rayleigh number and depth. Finally, we find that pseudo-steady linear and weakly nonlinear techniques described elsewhere are useful tools for predicting the behavior of instabilities beyond their formal range of validity, as good agreement is obtained with the simulations.

Feedback control of unstable cellular solidification fronts.

We present a feedback control scheme to stabilize unstable cellular patterns during the directional solidification of a binary alloy. The scheme is based on local heating of cell tips which protrude ahead of the mean position of all tips in the array. The feasibility of this scheme is demonstrated using phase-field simulations and, experimentally, using a real-time image processing algorithm, to track cell tips, coupled with a movable laser spot array device to heat the tips locally. We demonstrate, both numerically and experimentally, that spacings well below the threshold for a period-doubling instability can be stabilized. As predicted by the numerical calculations, cellular arrays become stable with uniform spacing through the feedback control which is maintained with minimal heating.

Chemoconvection patterns in the methylene-blue-glucose system: weakly nonlinear analysis.

The oxidation of solutions of glucose with methylene-blue as a catalyst in basic media can induce hydrodynamic overturning instabilities, termed chemoconvection in recognition of their similarity to convective instabilities. The phenomenon is due to gluconic acid, the marginally dense product of the reaction, which gradually builds an unstable density profile. Experiments indicate that dominant pattern wavenumbers initially increase before gradually decreasing or can even oscillate for long times. Here, we perform a weakly nonlinear analysis for an established model of the system with simple kinetics, and show that the resulting amplitude equation is analogous to that obtained in convection with insulating walls. We show that the amplitude description predicts that dominant pattern wavenumbers should decrease in the long term, but does not reproduce the aforementioned increasing wavenumber behavior in the initial stages of pattern development. We hypothesize that this is due to horizontally homogeneous steady states not being attained before pattern onset. We show that the behavior can be explained using a combination of pseudo-steady-state linear and steady-state weakly nonlinear theories. The results obtained are in qualitative agreement with the analysis of experiments.